Magnetic field induced directional localization in a 2D rectangular lattice 
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We study the effect of a perpendicular uniform magnetic field on the dissipative conductivity of a 
rectangular lattice with anisotropic hopping, tx ^ ty. We show that the magnetic field may enhance 
dramatically the directional anisotropy in the conductivity. The effect is a measurable physical 
realization of Aubry's duality in Harper systems. 
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Localization in aperiodic systems has been at the cen- 
ter of attention for decades in condensed matter physics. 
Most work concerned disordered systems. Twenty years 
ago, however, Aubry et al. [Q predicted that a ID tight 
binding hamiltonian with a quasiperiodic potential ex- 
hibits a metal-insulator transition as the amplitude of the 
potential becomes larger than a critical value. The proof 
rests on a duality property that allows mapping low into 
large coupling constants, while corresponding extended 
wave functions are transformed into states that are lo- 
calized. The tight binding model used by Aubry leads 
to the almost Mathieu equation |g,^ , a special case of a 
more general class of quasiperiodic systems for which the 
duality applies. It also happens that the almost Mathieu 
equation arises in the study of the dynamics of electrons 
in 2D in the presence of a rectangular lattice and a per- 
pendicular uniform magnetic field (see H for the case of 
a square lattice). In this case Aubry's duality may be 
interpreted as a rotation by 7r/2 of the lattice, an oper- 
ation that can be performed easily in a real sample and 
thus lends itself to experimental test. The aim of this 
letter is to show that due to Aubry's duality, turning on 
a magnetic field may produce a dramatic enhancement 
of the anisotropy already present in the conductivity of 
a rectangular potential. We denote by tx,ty the hop- 
ping amplitudes along the x and y axis (with tx < ty), 
and by <Jxx,<yyy the corresponding longitudinal conduc- 
tivities. At zero magnetic field the Drudc formula away 
from the parabolic edges of a tight binding band yields 
[pxxjoyy)^ K, const[txlty) ioi Small tx/ty. With mag- 
netic field, however, we obtain within the relaxation time 
approximation (RTA) 
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for irrational flux per plaquette. Here n(/i) is the density 
of states per unit cell at the Fermi energy ^, t is the 
dominant scattering time and 7 is a constant of 0(1). In 
deriving this result we assume that the temperature is 
low enough so that kT < Tl/t, and that the Fermi level is 
not too close from the edge of a gap larger than or equal 



to 0{kBT,h/T). It may fail as well in a regime where 
Mott's variable range hopping dominates. The enhanced 
asymmetry exhibited by Eq. (yj) for a large relaxation 
time is physically understandable in terms of Aubry's 
duality: the electronic eigenstates are extended along the 
easy direction y leading to a metallic-like behaviour for 
Gyy, whereas they are localized in the direction x, leading 
to an insulating-like behaviour for axx- We shall argue 
that the predicted enhancement should be observable in 
a superlattice of quantum dots. 

Consider a tight-binding hamiltonian in the x-y plane. 
For convenience we choose the gauge Ax = i?(0, x, 0) in 
which the coordinate y becomes cyclic, permitting plane 
wave solutions along this latter axis. The wave function 
along X must then obey the almost Mathieu difference 
equation |5| 



2ty cos{2'Kam — kyb)'ip{ma) 
tx [ip{{m + l)a) + ilj{{m ~ l)a)] = E%}]{ma) 
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Here the field variable a — eBab/hc gives the number of 
flux quanta traversing the unit cell, E is the energy, a 
and b are the lattice constants in the x and y direction, 
respectively, ky is the wave number of the free running 
plane wave along the y axis, and m is an integer labelling 
the lattice sites. The conventional Harper model y] is 
obtained by making the lattice square, with a = b and 
tx — ty. We are interested in the asymmetric case, usu- 
ally arising from the unit cell being rectangular, although 
a square array of elliptical quantum dots, for example, 
would also provide the required asymmetry. 

Inspection of Eq. ^ shows that in the limit tx/ty ^ 1, 
the solutions are plane waves, slightly modulated by the 
quasiperiodic potential. In the other extreme tx/ty ^ 1 
however, the solutions are localized features a distance qa 
apart if a = p/q is rational, or a single localized feature 
if this parameter is irrational 0. For tx/ty finite one can 
show that under Fourier transformation Eq. (0) formally 
turns into itself, with the roles of tx and ty exchanged. 
An extended state obtained for tx/ty ^ 1 is thus re- 
placed by a localized state in Fourier transformed space. 
This property is known as Aubry's duality [Q. Another 



way of obtaining an exchange of roles of t^ and ty is to 
change the gauge. To see this assume that t^/ty <^ 1 
so that the states given by Eq. (1) are locahzed along 
the X axis. If in the original problem one uses the gauge 
Ay = B{—y, 0,0) instead of A^ — B{0,x,0), the result- 
ing equation is formally identical to Eq. (1), only that 
now the wave function describes the dynamics along y, 
and tx,ty exchange places. Because of this latter fact and 
our assumption about the relative size of these parame- 
ters the new version of Eq. (1) gives now extended states 
that run along the y direction. Thus, while one gauge 
yields states that are localized along the hard hopping 
direction, the other gauge yields extended states along 
the easy hopping direction. This means that Aubry's 
duality is manifested in a single sample, its two principal 
axes playing the role of the dual states. As we show be- 
low, the anisotropy in the conductivity may reveal this 
effect in a dramatic way. 

The spectrum of Eq. (||) for irrational values of the 
field parameter is a Cantor set with a hierarchy of gaps 
that become extremely small [HJSJjg]. Fig. |l| shows the 
spectrum for different values of tx/ty. Note that as this 
ratio decreases from 0.8 to 0.2 the spectral gaps become 
more and more invisible due to their diminishing size. 
Each gap is uniquely labelled by an integer s ^ taking 
values between — (g — l)/2 and (g — l)/2 if a = p/q, and 
all values if a is irrational. It has been shown that the 
gap width behaves like 0] 
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for small tx/ty. For a rational a = p/q, the spectrum 
has exactly q subbands that do not overlap so that up 
to g — 1 gaps may appear. These are actually all open 
except for the one at £' = for q even ^^. For irra- 
tional a the spectrum can be well approximated by the 
rational approximants Pn/ Qn l9| obtained from truncat- 
ing the continuous fraction expansion of a at its n-th 
step ||l^ . The gap labels are then stable through this 
approximation 
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FIG. 1. Spectra of the Harper problem on a rectangular 
lattice for t^/ty = 0.8,0.6,0.4 and 0.2 obtained for rational 
values of a = p/q with q < 37. 

In real samples the presence of scattering limits the ex- 
perimental access to such a fine structure [Q. There 
is thermal broadening of size ksT, and also an energy 
width h/r associated with other sources of scattering, 
with T a characteristic relaxation time. All measurable 
quantities will be rather insensitive to gaps smaller than 
S — max (fc^r, ?i/T). Moreover, only energies within an 
interval of size S from the Fermi energy /i will matter for 
the electronic transport. Therefore, if a is irrational, it 
is sufficient to replace it by its best rational approximant 
p/q such that the gaps closest to fi have width not smaller 
than S. The error introduced in transport properties by 
this substitution is of the order of the hopping probabil- 
ity between sites at distance qa apart, which for tx small 
is bounded to be also very small. In our calculation of 
the conductivity we shall adopt this simplifying criterion. 
More precisely, for each integer n let A^") be the min- 
imal width of the closest gaps to /x corresponding to the 
rational approximant Pn/qn, and let s„ be the corre- 
sponding gap label. Then one chooses the largest value 
A^ of 71 such that 



A(JV+i) < g 



:(fcBT,^)<AW, 



(4) 



thus fixing the values oi q = q^ and s = sn- In par- 
ticular, as the temperature T decreases, one expects r 
to increase, so that the value of N increases as well and 
with it, those of qN and sn- 

A convenient form of Eq. (||) for the rational case is 

2ty cos(27rm- - kyb)<j)i{m) + tx [e'''^''(f>e{m + 1)+ 

e-^'^'-^f (m - 1)] = Ee{kx,ky)Mm) , (5) 

where m is an integer, (jyeim) is periodic of period q, and 
kx is a phase. For each point k — {kx, ky) in phase space, 
with < kxO,, kyb < 27r, there are q eigenvalues which, as 
k covers its range, make up the q subbands labelled by 
the integers ^ = 1,2, ...,g. All eigenvalues are contained 
in the energy interval \Ee\ < 2{tx + ty) ~ W/2, with W 
the width of the original zero-field band. We assume the 
Fermi energy /i to lie somewhere within this range. 

In the infinite volume limit and in the relaxation time 
approximation, the longitudinal conductivity for our q 
subbands system, is given by Kubo's formula |I3] 
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47r2 \/T~i{Ei~Et') /fi 
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where Eg and E^i depend upon the phase space coor- 
dinates {kx,ky) and the integrals are taken over all of 
phase space. Following condition (Q) , in Eq. (g) we have 
replaced the Fermi distribution by a step function. The 
first term in this equation is the intraband term, while 
the second includes all interband contributions. In the 
absence of a magnetic field there is just one band in our 
model and only the first term is relevant. In what follows 
we shall assume t^ <ty 

A careful analysis of Eq. (ph shows that the conductiv- 
ity depends strongly on the position of the Fermi energy 
with respect to subband edges. We distinguish two cases. 
1) n far from the subband edges: assume first that 
the Fermi energy is at a distance larger than O (A^/f^) 
from the nearest gap. The intraband term can be com- 
puted solving the eigenvalue equation 



Pq{E)^2{~Y+^ {tl COS qat^ + tl COS qbky 



(7) 



for the appropriate subband £. In this expression Pq is the 
Chambers polynomial associated with Eq. (|^) ||l^,01. Due 
to condition (0) , the interband contribution in Eq. (O) can 
be expanded in 1/r since Ti/t < \Ei — Ee>\. Keeping just 
the first terms in this expansion, one gets to lowest order 
in A = t^/ty 
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Here, as mentioned earlier, n(^) is the density of states 
at the Fermi energy and 7 is a numerical factor of or- 
der 1, while r = a/b and <^ = Pq{^)/2t'^. The first term 
in these expressions is the Drude conductivity coming 
from the intraband transitions in Eq. (j^) . That the ratio 
{'^xxl(yyv)Y)^^^c is °f 0(A^'') is apparent from the energy 
derivative in Eq. (g) and the fact that the Chambers 
polynomial depends on the phases through the constant 
term in the right hand side of Eq. (0), only. Also, that 
the lowest order contribution to the interband transitions 
is O(A^) follows from the fact that in a perturbative ex- 
pansion in terms of A, the zcroth order term in <j)i does 
not depend on the phases. 

Assuming g^r > 1 and having in mind that at suffi- 
ciently low temperatures and within the rational approx- 
imation ansatz \sn\ < (gAT — 1)/2 and Ti/rij, ~ A'^^^/ij, ~ 
X\sn\ > ^(9iv-i)/2^ Qjjg finds the ratio between the Drude 
and interband contributions to the conductivity along the 
X direction to be negligible. Ignoring the first term in dq) 
we then have 



h ■Krn{n) 
1 I 7, — ^ 






(10) 



This important result states that as long as ^ is not too 
close to a subband edge of the proper rational approx- 
imant, the conductivity ratio vanishes quadratically in 



the inverse relaxation time for small A. Applying Eq. (6) 
to the original tight binding band to obtain the corre- 
sponding conductivity ratio in the absence of a magnetic 
field near the band center one then arrives at the form 
given by Eq. (1). 

2) /i near a subband edge: when the Fermi energy /i 
is very near a subband edge the situation changes signifi- 
cantly. The density of states has a logarithmic singularity 
in that neighborhood, and n(/^) oc (iy/ix)'"^^ in a region 
of order A^/iy from the edges. As a result, Eqs. (4)- 
(5) are no longer correct. The expressions for the Drude 
contribution are valid in the region < |^| < 1 — A'" 
only. In the range 1 - A«" < |^| < 1 -|- A«« or if the 
Fermi energy lies in a gap, a^x may become as large as 
(jyy . A detailed discussion of this regime will be published 
elsewhere ]1J]. These situations may not be relevant for 
experiments however, because they only occur in a very 
small region of the spectrum. Indeed, when t^ < ty, 
the probability of having // within O (A^/ij,) from the 
subband edge behaves like (tx/ty)"^. In addition, for qn 
large the total length of the subbands is given by 4:\ty—tx \ 
|l|, [l6|Jl7| . Thus, choosing the Fermi energy fi at random 
in the energy interval [—W/2, W/2] gives /i in one sub- 
band with relative probability 
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Thus, the probability of having /i in a gap or close to 
a gap edge is proportional to tx/ty which is small Q. 
Moreover, it is likely that very small fluctuations in the 
magnetic field will tend to wash away the effects of /i 
lying in the anomalous regions close to a subband edge 
or in a gap. 

In summary, we have shown that in the presence of a 
perpendicular magnetic field the field free conductivity 
asymmetry of a rectangular lattice may be dramatically 
enhanced in a pure sample. This is a physical realization 
of Aubry's duality that may be understood in the follow- 
ing way. The magnetic field affects the hopping of the 
tight binding electrons through the effective potential of 
period q times the appropriate lattice constant. The site 
energies are then no longer the same within and the elec- 
tron has to tunnel a distance equal to this period, through 
a potential mismatch that scales with A~^ or A, depend- 
ing on whether transport is along the x or y directions. 
For an irational field the tunneling can only be made 
possible through scattering events, which are less likely as 
the relaxation time diverges. Evidence for this behaviour 
has also been found in the spreading of a wave packet 
in the presence of a weak modulation potential in two 
dimensions. In the regime in which the single band ap- 
proximation holds, such spreading is entirely directional, 
occuring along the direction of largest hopping amplitude 
only pSl . Our results also shows that the suppression of 
the conductivity along one of the principal axes by an ex- 
ternal field is controlled by the ratio of the gap size to the 



zero-field band width (;i/r)(l/W) ~ (i^/tj,)!""!. With 
present day rectangular arrays of antidots this quantity 
is of order one tenth. Thus the effect should already be 
visible. It is likely that, as the miniaturization of meso- 
scopic technology progresses, this quantity can be made 
even smaller so that the magnetic field induced localiza- 
tion may be more easily observed. 
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